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Abstract 

Statistical fragment emission from excited nuclear systems is studied within the framework of a 
schematic Fermi-gas model combined with Weisskopf's detailed balance approach. The formalism 
considers thermal expansion of finite nuclear systems and pays special attention to the role of the 
diffuse surface region in the decay of hot equilibrated systems. It is found that with increasing 
excitation energy, effects of surface entropy lead to a systematic and significant reduction of effective 
emission barriers for fragments and, eventually, to the vanishing of these barriers. The formalism 
provides a natural explanation for the occurrence of negative nuclear heat capacities reported in 
the literature. It also accounts for the observed linearity of pseudo-Arrhenius plots of the logarithm 
of the fragment emission probability versus the inverse square-root of the excitation energy, but 
does not predict true Arrhenius behavior of these emission probabilities. 

PACS numbers: 21.65+f, 21.60.Ev, 25.70.Pq 
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I. INTRODUCTION 



Over the last decade, a considerable effort has been made to understand the phenomenon 
of production of multiple intermediate-mass fragments (IMF) in individual nuclear reaction 
events. Yet, in spite of this effort, understanding of this process is still far from being 
satisfactory. For example, a consensus is still lacking on even such a fundamental issue of 
whether the fragment production is dominantly a statistical or a dynamical process. To a 
certain extent, such lack of consensus is due to the fact that different models, promoting 
conflicting production scenarios, appear equally successful in describing select subsets of 
experimental observations. 

The present study is driven by the realization that statistical models that have been used 
in the interpretation of nuclear multifragmentation data Jl], |^, ^ |4|, [5[ fail to provide a coher- 
ent or satisfactory IMF production scenario. Among other things, they fail to overcome the 
fundamental difficulty of traditional equilibrium-statistical decay models |], [7]] in explain- 
ing how IMF production can effectively compete with nucleon and light-charged particle 
production, when the IMF emission barriers are many tens of MeV high, compared to ~ 6 
MeV for protons. It was the inability of these statistical decay models to explain the IMF 
production rates that gave rise to the attractive speculation that multiple IMF production 
could be a manifestation of the liquid-gas phase transition in finite nuclear matter. This is 
so, because one would expect a gas of nucleons enclosed in a container of sufficiently large 
size to equilibrate under certain circumstances into a mix of nucleons, light-charged parti- 
cles, and IMFs. This would be a process essentially analogous to condensation of a gas into 
liquid droplets. If this were the case, the abundance of IMFs could be large and, perhaps, 
comparable to that observed experimentally. 

For many years, speculations regarding the relevance of the liquid-gas transition for what 
is termed nuclear multifragmentation have driven studies of this latter phenomenon. So 
far, however, scrutiny of practical models featuring such a phase transition @, |], || 
has revealed significant gaps in the chain of reasoning leading to the "prediction" of an 
experimentally observable liquid-gas phase transition. For example, models based on the 
concept of an existence of a freeze-out configuration || |J utilize in their formalisms a 
"thought" container with ideally reflecting walls, the boundary of the freeze-out volume. 
Reflection of any outgoing flux of nucleons and light-charged particles (LCP) back into 
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the freeze-out volume, |8j is crucial for the models allowing particles to eventually form 
bigger clusters, the intermediate-mass fragments. However, since there is no known physical 
boundary defining such a freeze-out volume, or mechanism by which IMFs could pass the 
high saddle of emission barriers, the predicted IMF abundance must be considered an artifact 
of the modelling itself. 

Similarly severe are shortcomings of approaches utilizing percolation [5[] or Fisher's model, 
|| which have recently been applied || to nuclear multifragmentation. jllj Notably, these 
models do not even address the issue of the high Coulomb barrier inhibiting fragment emis- 
sion. The presence of a high barrier is obvious even in the vicinity of the critical point, where 
the expected matter density (in the range of 1/3 - 1/2 of ground-state matter density) is 
still high enough to trap any fragment formed inside the expanded nuclear system in very 
much the same way as it would trap fragments formed at ground state density. In the case of 
Fisher's model, largely overlooked remains the fact that, by design, the formalism is limited 
to events in which there is no more than one cluster present at any time, with the balance of 
the system being always in a state of a homogenous, monomeric gas of individual molecules. 
While such "either-zero-or-one-fragment" events may be essential in a condensation theory, 
such as represented by Fisher's model, they are of little, if any, interest in the context of 
nuclear multifragmentation. 

A somewhat different statistical approach to explaining high yields of IMFs observed in 
energetic nuclear reactions was adopted in the expanding-emitting source model (EESM). 
U This model considers dynamic, isentropic, self-similar expansion of a nuclear system to 
densities below the equilibrium density proper for the given total excitation energy, where 
it finds enhanced fragment emission probabilities. However, as pointed out recently, || 
the EESM formalism unjustifiably derives most of the (collective) energy necessary for a 
fragment to overcome the Coulomb emission barrier, not from a thermal heat bath, as 
prescribed by Weisskopf 's approach, but from a thermally isolated single, compressional 
degree of freedom. 

The above discussion emphasizes the need for an alternative scenario of statistical IMF 
emission. The present study demonstrates that a simple scenario, closely related to that 
known from fission studies, offers an explanation of how at moderately high excitation 
energies, IMF emission can compete effectively with nucleon evaporation. It represents 



an extension of ideas and formalism presented in a recent publication. [12] Central to 
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this formalism is the notion of a relatively high entropy associated with the diffuse nuclear 
surface region (as opposed to bulk matter). Note that it is also the increased surface entropy 
that enhances the fission probability of a heavy nucleus relative to nucleon evaporation. 
This latter enhancement is commonly described in terms of the ratio a//a n of level density 
parameters for saddle and ground-state configurations, respectively. 

The formalism adopted in the present study is described in detail in Section II. While this 
formalism is somewhat schematic, it is believed to capture the essential physics underlying 
the processes involved. One benefit of such a schematic treatment is that it provides a 
direct insight into the phenomena of interest, disregarding a multitude of secondary details 



demanded by a more rigorous approach. Results of calculations are presented in Section III 



while the conclusions are presented in Section IV 



II. THEORETICAL FORMALISM 

The present study assumes that an excited nuclear system expands in a self-similar fashion 
so as to reach a state of approximate thermodynamical equilibrium, where the entropy S is 
maximal for the given total excitation energy E^ ot , i.e., 

— dp — lE ^-°- (1) 

The functional dependence of entropy on E* ot and the bulk nuclear matter density p is 
evaluated using the Fermi-gas model relationship 



S = 2y aEt h = 2y a(E^ ot - E compr ), (2) 

where a is the level density parameter, E^ h is the thermal excitation energy and E compr is 
the collective compressional energy. The dependence of S on bulk nuclear matter density 
p arises in Eq. ^| through the dependence of both, the level density parameter a ( "little-a" ) 
and the compressional energy E compr on the matter density. 

The dependence of little-a on the nuclear matter density for infinite nuclear matter is 
given by the Fermi-gas model: 



Oo(f)-*, (3) 

Po 
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where a Q is the level density parameter for the nuclear matter at ground-state matter density 

Po- 

The above equation holds approximately also for finite nuclei if the expansion or com- 
pression of these nuclei is assumed to occur in a self-similar fashion. This is so because for 
finite nuclei, the little-a parameter consists ]TB|, |T^] of volume and surface terms, ay and a a , 



respectively, both of which are proportional to p 2 / 3 under the assumption of self-similar 
expansion: 

O 2 2 O 2 

a = a v + a a = A(—yza v + A*(—)~3a a} (4) 

Po Po 

where ay and a a are volume and surface coefficients, respectively, independent of bulk 
nuclear matter density. 

The term "self-similar expansion" is used here to describe a type of expansion in which 
any change in the matter density profile is reducible to a simple rescaling of the radial 
coordinate, such that: 



f p (r) = c 3 f (cr), (5) 

where / (r) is the ground-state density profile function and c is a scaling constant. 

The presence of a surface contribution to the level density parameter is of crucial impor- 
tance in the present study as it describes the part S a of entropy S of the system associated 
with the diffuse surface domain and is seen to have significant effect on the fragment emission 
probability. One has 

S = Sy + S a = a v T + a a T, (6) 

where Sy is the entropy of the bulk matter and T is the system temperature. 

The compressional energy in Eq. ^| is approximated in the present study following the 
schematic prescription proposed in the expanding emitting source model EESM |], i.e., 

^compr Q>(1 ) ' (^0 

Po 

where e compr and e& are the compressional and the ground-state binding energy per nucleon 
of the system, respectively. Note that Eq. [1] ensures that the compressional energy varies 
parabolically with p, from zero at ground-state density (p G ) to e b at zero density. 
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Equations [T|-|3] and [7] allow one to obtain an analytical expression for the equilibrium 

density, p eq / p Q , of nuclear matter as a function of the excitation energy per nucleon, e\ ot = 
Et ot /A, where A is the mass number of the system: 



^ = 7(1 + 4/9-8^). (8) 

Po 4 V e b 

Equation [8] reflects the fact that for e\ ot < e b , the system is bound as far as the self- 
similar expansion mode is concerned, and features a single maximum in entropy in the 
range of matter densities 1/2 < p eq j ' p Q < 1. For e\ ot > e b , the system is essentially unbound 
with respect to the self-similar expansion mode, with the entropy diverging as the matter 
density p tends to zero. However, for (9/8)e;, > e* tot > e b , the system still has a local 
metastable maximum at a finite density given by Eq. || i.e., in the range of matter densities 
1/4 < Peq/ Po < 1/2. The latter metastable (with respect to self-similar expansion mode) 
maximum in entropy is separated from the divergence at zero density by a minimum at 
P 'saddle, where 



P^dMe = h 1 _ /g _§£*£*). ( 9 ) 

Po 4 V e b 

Here, E b = Ae^- 

Further, one can obtain the caloric equation of state of an equilibrated system (i.e. at 
constant, zero external pressure) by using the Fermi-gas model relationship between the 
system temperature T and thermal excitation energy E^ h 



T = \P± = {^a-\\Et ot -E b {l- P -^y. (10) 
V a Po V Po 

The probability p of emitting a fragment from an equilibrated excited system as defined 
above can be evaluated using the Weisskopf formalism: [TT| 



p oc e = e saddle eq , (11) 

where S sa ddie and S eq are saddle-point and equilibrium-state entropies, respectively. The 
latter two entropies can be calculated using Eq. 0: 



S eq = 2 a A [E* tot -E b {l-^)\ (12) 

V Po 
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Ssaddle — S res + Sfrag — 2y (a res + 0,f rag )E* addle . (13) 

In Eqs. [12| and a a, a res and a frag are the level density parameters of the system 
at equilibrium, the residue, and the fragment, respectively, while E* addle is the thermal 
excitation energy of the system in the saddle-point configuration. The latter quantity is 
calculated as 

E *sa H ddle = E tot ~ E b (l- —) 2 - V sadd l e , (14) 

Po 

where V sa ddie is the (collective) saddle-point energy. 

Formally, the emission probability p can be also expressed in terms of an effective emission 
barrier B e ff, 

p oc e T , (15) 

when one sets 



B eff = -TAS. (16) 

A selection of results of calculations performed using the above formalism is presented in 
Section |TJ below. 



III. RESULTS OF MODEL CALCULATIONS 



Results of the calculations performed in the framework of the formalism presented in 
Section II are presented in Figs. 1-7. In these calculations, values of ay =1/14.6 MeV -1 
and a CT =4/14.6 MeV" 1 have been assumed for the coefficients ay and a a , as suggested in 
literature []I3"| . Further, e =8 MeV was assumed [|I|] for the ground-state binding energy 
per nucleon, while the saddle-point collective energy was approximated by the Coulomb 
energy of the residue and fragment represented by two touching spheres of radius parameter 
rcoui = l-3(p/po) -1 / 3 fm. The calculations were made for excited 197 Au nuclei. 

Fig. Q illustrates the dependence of the equilibrium density, p eq , of bulk matter on the 
excitation energy per nucleon (solid curve). As seen in this figure, for the range of excitation 
energies readily accessible in experiments, the bulk matter density in a state of maximum 
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FIG. 1: Equilibrium density of bulk nuclear matter (solid curve) and the net gain in entropy 
realized by an 197 Au system as a result of the relaxation of the self-similar expansion mode (dashed 
curve), plotted as functions of total excitation energy per nucleon. 

entropy differs substantially from that of the nuclear ground state. This affects both, the 
caloric equation of state and the fragment emission probability, two entities of considerable 
interest . 

The dashed curve seen in Fig. p] illustrates the net gain in entropy resulting from the 
relaxation of the self-similar expansion mode in an excited 197 Au system. Large gains in 
entropy associated with the relaxation of this mode emphasize the importance of this mode 
for a statistical description of excited nuclear systems, notably for models based on the 
concept of microcanonical || or pseudo-microcanonical Q equilibrium. 

The caloric curve calculated for states of maximum entropy, i.e., for the states of equi- 
librium density p eq is depicted in Fig. [| Not surprisingly, this curve shows considerable 
deviation from the simple Fermi-gas form of T oc \fEfot- Note that in experiments, such as 
the recently reported ISIS experiment, (jDJ it is E^ ot and not the purely thermal contribu- 
tion to it, E* h , that is in fact measured. This is so, because the static compressional energy 
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FIG. 2: Caloric curve for the 197 Au system. 

Ecompr is experimentally undistinguishable from thermal excitation E^ h . 

Most notably, the caloric curve seen in Fig. |2| features a maximum temperature of approx- 
imately T max =Q MeV, an indication that for higher temperatures the system is inherently 
unstable and does not find an equilibrium density. In other words, under the assumption that 
its matter distribution is homogeneous, a nuclear system placed in a heat bath of T > 6MeV 
would expand indefinitely, for which it would derive increasing amounts of energy from the 
heat bath. In a more realistic case, which is beyond the present consideration, the sys- 
tem would likely decay into a "gas" of fragments and free nucleons that would continue its 
indefinite expansion. 

The particular form of the caloric curve seen in Fig. |^ can, perhaps, be better understood 
when inspecting the dependence of the free energy F on matter density p and temperature 
T, which is is illustrated in Fig. || The calculations show that, at temperatures below 
T « QMeV, the free energy as a function of p features two equilibrium points, a local 
minimum and local maximum, where the latter reflects an unstable equilibrium These two 
points predicted for any temperature T of less than approximately 6 MeV, correspond to 
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FIG. 3: Free energy for the Au system as a function of nuclear matter density and temperature. 

different excitation energies per nucleon, e* tot . Note that at both these extremal points, the 
pressure of the system at its surface is zero. This is so, because the pressure is proportional 
to the partial derivative of the free energy with respect to density at constant entropy, and 
both, the free energy and the entropy are stationary with respect to the matter density at 
the extremal points in question. At T ~ 6MeV, the minimum and maximum in F as a 
function of p merge into an inflection point and, then, at even higher temperatures, the 
free energy features only a monotonic decrease with decreasing p. The inflection point in F 
corresponds to maximum in T as a function of E^ ot . Obviously, a monotonic decrease in free 
energy with decreasing matter density means that no equilibrium density exists for a given 
temperature. 

The role of the surface entropy in fragment emission is clear from Fig. |] illustrating the 
dependence of effective barriers (see Eq. [16]) for the emission of 12 C and 16 O from exited 197 Au 
nuclei. Solid lines in this figure illustrate the results obtained when the surface entropy is 
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FIG. 4: Effective barriers for the emission of C and O fragments from equilibrated excited 



taken into consideration via inclusion of a a in the expression for the level density parameter 
a (see Eq. |3|), while dashed lines represent results of calculations in which the surface entropy 
effects were neglected, i.e., in which a was assumed simply proportional to the mass number. 

It is important to note that with the reduction in the effective barrier being dominantly 
due to the surface entropy effects, no analogous reduction is expected for the emission of 
nucleons or light charged particles. In the latter case, the effective barriers are expected to 
show trends similar to those shown by the dashed lines in Fig. f|. 

The effects of the surface entropy on the relative emission rates of various IMFs are 
illustrated in Fig. |5|, displaying the reduced probabilities (taken as bare Boltzmann factors) 
for the emission of IMFs of different atomic numbers from an equilibrated 197 Au system 
as functions of the total excitation energy per nucleon. The spectra are normalized to the 
probability for emitting a proton, bound to the system by Q p =8 MeV and experiencing a 
Coulomb barrier of 4 MeV (corresponding to a Coulomb radius parameter of rc=1.3 fm). As 
seen in this figure, already at excitation energies per nucleon of the order of 3 MeV/nucleon, 



197 Au systems as functions of excitation energy per nucleon. 
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FIG. 5: Distributions of relative emission probabilities of various IMFs from an excited, equili- 
brated 197 Au system as functions of the total excitation energy (solid curves). The dashed line 
represents the boundary of the domain of dynamical instability of the system. (See text) 

fragment emission begins to compete effectively with proton emission and, then, at higher 
excitation energies, fragment emission becomes the dominant decay mode. Note, that the 
highest probabilities seen in Fig. [|, above the dashed line, reflect dynamical instability with 



respect to the fragment emission mode. In fact, Weisskopf's approach is inapplicable in 
the domain above the dashed line, as in this domain, the saddle-point entropy exceeds that 
for the "equilibrium" density of a self-similarly expanded system. In the calculations, it was 
assumed that the mass number of the fragment is twice the atomic number, i.e., Af r = 2Zf r . 

Since the magnitude of the effective barrier depends on the excitation energy and, thus, 
on temperature T, and since T is a non-monotonic function of E^ ot , one would expect the 
fragment emission probabilities p to deviate significantly from a simple Arrhenius law with 
an exponential functional dependence on the inverse temperature 1/T. This expectation 
is confirmed by Fig. 0, where logarithmic plots of p versus 1/T are seen to feature clear 
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FIG. 6: Arrhenius plots for emission of C and fragments from excited Au systems. 

deviations from linearity, including a prominent "back-bending". The latter back-bending 
is obviously expected in view of the maximum in the caloric curve seen in Fig. Such a 
behavior is also suggested by the general deviation of the caloric equation of state from the 
one for a low-temperature Fermi gas. 

On the other hand, pseudo- Arrhenius plots depicted in Fig. |7|, where the logarithm of the 
emission probability, ln(p), is plotted versus the inverse square root of the total excitation 
energy, 1 / y/Et ot , , are to a good approximation straight lines. This observation comes as a 
surprise, as it cannot readily be expected based on the details of the theoretical formalism 
employed. Therefore, no simple explanation for such a linearity can be offered at this time, 
the very rationale behind the construction of such pseudo-Arrhenius plots being numerous 
experimental findings. ||, [16], |I7|, [T!| [L9| While it may be purely fortuitous, the linear char- 
acter of pseudo-Arrhenius plots predicted by the present formalism seems to be in agreement 
with experimental observations. 
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FIG. 7: Pseudo-Arrhenius plots for emission of C and fragments from excited Au systems. 
IV. CONCLUSIONS 

A formalism has been developed to describe quantitatively, albeit in a schematic fashion, 
a scenario of purely statistical emission of massive fragments from finite equilibrated systems. 
The formalism recognizes the importance of thermal expansion of hot matter and considers 
the stability of such systems at equilibrium nuclear densities. The important role of surface 
entropy consists in an effective "softening" of the nuclear surface, resulting in enhanced 
fragment emission probabilities. The decay scenario underlying this formalism is that of 
a thermally expanded system with developed thermal fluctuations of the diffuse surface. 
Within such a scenario, fragment emission is naturally viewed as analogous to fission. 

The calculations performed within the framework of the proposed formalism show that, 
with increasing excitation energy and mostly due to effects of increased surface entropy, 
the effective barriers for fragment emission decrease and, eventually, vanish. As no equiva- 
lent decrease is expected for nucleons and light charged particles, one expects that at high 
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excitation energies, fragment emission competes effectively with nucleon and LCP emission. 

The formalism presented here also predicts an occurrence of negative heat capacity in 
highly excited nuclear systems. It also predicts linear pseudo-Arrhenius plots, where the 
logarithm of the fragment emission probability, ln(p), is plotted versus the inverse square- 
root of the excitation energy, \j\fE*. This latter prediction is found in agreement with 
numerous experimental observations. [|ll], [T7|, |TB], However, as also expected, the formal- 
ism predicts non-linear true Arrhenius plots, where the independent variable is the inverse 
temperature, 1/T, rather than 1/yE*. This is so, because the effective emission barriers 
depend rather strongly on excitation energy (and temperature) and, hence, fail to satisfy 
one of the fundamental prerequisites of the Arrhenius law. 

The simplicity of the presented formalism should allow one to incorporate it into existing 
statistical decay codes. For purposes where the ultimate goal is to obtain an agreement 
with experimental data, it would be desirable, however, to improve the formalism itself. 
For example, one should use more realistic (and lower) saddle-point energies, instead of the 
schematic Coulomb energies of touching spheres. One should also be mindful of differences in 
equilibration times for nucleon and LCP emission modes, on the one hand, and for different 
IMF emission and fission modes, on the other hand. It is natural to expect that modes 
requiring major rearrangement of mass require longer times to equilibrate than modes that 
do not require such a rearrangement. Consequently, one would expect all IMF emission 
modes to attain equilibrium later than nucleon emission modes, and fission modes to reach 
equilibrium even later. As a result, early in the decay process, while the IMF degree of 
freedom is approaching equilibrium, one would expect to see "pre-IMF nucleons" , in much 



the same way as one sees pre-fission neutrons. p(| [HJ] Similarly, one would expect to see 
pre-fission IMF emission in situations when for the same systems at equilibrium, fission 
would dominate. Importantly, the latter pre-fission IMF emission modes may succeed in 
removing enough charge and excitation energy from the system so as to pre-empt fission 
altogether - a phenomenon that has recently been observed experimentally. 
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